What is the Thouless Energy for Ballistic Systems? 
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The Thouless energy, E c characterizes numerous quantities associated with sensitivity to bound- 
ary conditions in diffusive mesoscopic conductors. What happens to these quantities if the disorder 
strength is decreased and a transition to the ballistic regime takes place? In the present analysis 
we refute the intuitively plausible assumption that E c loses its meaning as an inverse diffusion time 
through the system at hand, and generally disorder independent scales take over. Instead we find 
that a variety of (thermodynamic) observables are still characterized by the Thouless energy. 



Much of the physics of diffusive disordered conductors 
is governed by an energy scale E c , known as the Thouless 
energy [1]. By definition, E c = TiD / L 2 (D is the diffusion 
constant and L the system size) coincides with the inverse 
of the diffusion time through the system, the latter be- 
ing a purely classical quantity. Nevertheless, there are 
numerous energy scales, related to quantum-mechanical 
aspects of disordered conductors, which equal (up to nu- 
merical factors) E c as well: 

E\\ The disorder averaged curvature (at $ = 0) of the 
total energy of an AB-system containing a constant 
(i.e. flux independent) number of particles. This 
quantity is proportional to the flux derivative of the 
(canonical) persistent current (I(<&)) d L near zero 
flux. Here (■■■) dL denotes impurity and sample 
size averaging. 

E 2 : The energy scale governing the suppres- 
sion of the harmonics content I m := 
2 ^o" 1 So" d$sin(27rra$/^ )/($) of the total per- 
sistent current (I^))^^ {4>o = h/e). According to 
Ref. [2], 

< 7 ™) d , £ = -4A/0Q exp (-2m V / 7A^) , (1) 

where 7 is an inelastic broadening [3], and A the 
average level spacing. 

E3 : The typical current carried by a single electron level 
(i.e. the derivative of the latter with respect to $) 

averaged over (« 2 ($)) ~ <t>o AE C . Here (. . .) 

\ / d,L 

refers to flux averaging. 

There are further quantities that are not manifestly 
boundary condition sensitive but nevertheless equal E c 
in diffusive systems: 

E4: the quantity g(h/e 2 )A where g is the average two 
lead conductance of diffusive conductors [1]. 



E$: The energy scale limiting the applicability of ran- 
dom matrix theory (RMT) to the description of 
spectral correlations. 

The above list is not exhaustive. Yet, all other occur- 
rences of E c known to us are closely related to at least 
one of the examples E\ , . . . , E$ . The realization of the 
important role played by the Thouless energy can be re- 
garded as one of the major steps towards a comprehensive 
theory of disordered systems. 

The above statements apply to diffusive metallic sys- 
tems, i.e. systems larger than the elastic mean free path 
I. Upon reducing the amount of disorder (or decreas- 
ing L) a transition to the ballistic regime takes place. 
This, however, does not imply that disorder becomes im- 
material. A simple heuristic argument supporting the 
opposite may be given as follows: Most of the above 
mentioned quantities are closely related to the system's 
quantum-mechanical energy spectra. In a generic d- 
dimensional ballistic system I > L > l(ppLy~ d (pp is 
the Fermi momentum) [4] the typical value of the dis- 
order scattering matrix element exceeds the mean level 
spacing by a large factor (ppL) d ~ 1 L/l which means that 
the distribution of the energy eigenvalues differs dras- 
tically from the clean case. This implies the existence 
of a distinct regime (rather than a crossover point), in- 
termediate between the diffusive and the clean limits 
[5]. The characteristic energy scales in this regime sat- 
isfy A < t -1 < tj 1 < E c <C Ep, where r is the 
elastic scattering time, tf = vpL^ 1 the ballistic time 
of flight, vp the Fermi velocity, Ep the Fermi energy, 
E c = vpl / (d- L 2 ) is defined formally as before and ft has 
been set to unity. The above hierarchy is to be compared 
with A < E c < t _1 <C Ep characteristic for the diffusive 
regime. 

The existence of strong spectral correlations suggests 
that at least some of the scales Ei, . . . , E 5 continue to 
be disorder sensitive even in the ballistic regime. On the 
other hand, E c loses its meaning as an inverse transport 
time through the system, in which respect it is replaced 
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by the - disorder independent - inverse time of flight tf . It 
is thus imperative to ask whether it is possible to find an 
energy scale in the ballistic regime playing role analogous 
and as crucial as E c in the diffusive case and whether this 
scale is disorder dependent. The results of our analysis, 
reported here, show that for boundary condition (or mag- 
netic flux) sensitive quantities the Thouless energy, E c , 
is still the relevant scale. 

To be specific, we consider an ensemble of two- 
dimensional AB cylinders of height and circumference L 
which differ in their impurity configuration and size. The 
relevance of sample size averaging will become evident 
below. Apart from the average conductance, E4, all the 
above listed quantities can be calculated by means of the 
correlation function 

K e ™(e, 4> 1 ,<f> 2 ) := \{v(E F + e, fa)v(E P , fa)) L d - 1 

characterizing fluctuations in the energy spectra of the 
individual systems. Here <f> := 2tt^/4o, v{e,<ft) is the 
flux dependent densitiy of states and vq := (v(e, 4)) l d is 
practically constant. In order to evaluate this expression 
in the particular case of a ballistic ensemble, it is essen- 
tial to note two features of the combined disorder/sample 
size average: (i) Size and disorder average are essentially 
statistically independent, {. . .) L d ~ {{. . ■) d ) L [4] and (ii) 
size averaging is equivalent to averaging over the Fermi- 
energy (see below). For our purposes it suffices to con- 
sider an average over a range E > Ai ~ (ppL) d ~ 1 A cor- 
responding to size fluctuations on scales ~ 1/pi?- From 
(i) and (ii) it follows immediately that [11] 

K^(e,4>i,4>2) = {K(e,4>u4>2)) Eo + K'^ufa), (2) 
where 

K'(e,4>u4>2) := X{{v{E F + *,fa))MEF,<h))d) e* ~ ^ 

K(e, ^,4)2) ■■= \ ({v{E F + e,4> 1 )v(E F ,4>2))d ~ 

-{v(E F + e,4>i))MEF,4>2)) d ), (3) 

In the diffusive regime, (v(e,4>))d — ^0 =^ — 
implying that the structure of the spectrum is solely de- 
scribed disorder induced correlation function K. On the 
other hand, it is clear from the defining equations that 
K — > in the clean limit, in which case the function 
K' , corresponding to correlations in the clean spectrum, 
gains importance. This motivates the above definition of 
the ensemble average: The structure of the clean spec- 
trum depends sensitively on various parameters, such as 
the system size or the position of the energy window un- 
der consideration. Any realistic modelling of an ensemble 
average should thus involve a simultaneous average over 
both, disorder and these latter characteristics. In an en- 
semble of diffusive systems the fluctuations induced by 
disorder are strong enough to mask any underlying clean 



structure and averaging over other parameters is inessen- 
tial. In the course of a crossover to the clean limit the 
structure of the clean spectrum becomes visible and the 
function K' begins to play a dominant role. The ballis- 
tic case represents a transition regime in the sense that 
the functions K and K' are simultaneously relevant. We 
next discuss both contributions separately. 

Within the framework of diagrammatic perturbation 
theory, the disorder induced correlation function K [4] 
can be written as the sum of a ladder type diagram 
and its time reversed analogue [3], here referred to 
as the diffuson and the Cooperon contribution, respec- 
tively. In the time reversal invariant (field free) case, 
both diagrams are equal and one arrives at if(e, 0, 0) = 
A 2 /(27r) 2 Re d 2 2£ 9 m(l - X(q,e)) + ..., where J2 g de- 
notes a summation over a two-dimensional vector q = 
(2im x jL x , -niiy/Ly), n x G Z,n y £ J\f of quantized 

modes, X(q,e) = [(1 + 7T + ier) 2 + (\q\l) 2 ] 2 , 7 ~ A 
is an inelastic broadening, ln(l + z) := ln(l + z) — z and 
the dots indicate small and field insensitive corrections 
which may be discarded in the present discussion. 

Upon the inclusion of magnetic fluxes, the de- 
generacy between the diffuson and the Cooperon 
contribution is lifted and we obtain K(e,4i,42) = 
A 2 /(2 7 r) 2 Re d 2 e F(e,4 u 4 2 ), where F(e,4> 1 ,4 2 ) = 

^ E„ [1 " A (q + L-^+e-,6)] + (4+ +>4>-)), 
where 4>± := 4i =t 42, 4 is the unit vector in re- 
direction and 4+ (4-) corresponds to the Cooperon 
(diffuson) contribution. For practical convenience we 
next introduce the function (SN 2 (4)) ens describing the 
flux sensitive contribution to the fluctuation of the to- 
tal number of occupied levels. This quantity relates to 
the above introduced function K ens via (<5Af 2 (</>)) ens = 

A- 2 /^ de x de 2 [K^ s {e 1 - e 2 ,4,4) - Jf ens ( £l -e 2 ,0,0)]. 
In analogy to Eq. (2) we further define (SN 2 (4)) ens =■ 
(8N 2 (4)) + {&N 2 (4))\ where the primed (unprimed) 
term corresponds to the K' {{K) Eq ) contribution to 
K ens . Due to the fact that K has been expressed in 
terms of a twofold derivative, the energy integration nec- 
cessary to obtain (SN 2 (4)) can trivially be performed 
and we arrive at {SN 2 (4)) = F(E F ,4,4) - F(Q,4,4) - 
F(E F ,0, 0) + f(0,0,0). An approximate harmonic ex- 
pansion of (SN 2 (4)) reads 

{SN 2 (4)) -^jt {-^ e ~ 2m ^~ c - Ae- 2mB '^\ ■ 

m=l ^ ' 

■ [cos (2m<j>) - 1] , (4) 

where A, B are constants of 0(1) (Within our numerical 
accuracy we obtain A ~ 0.92, B ~ 1.21). 

We next discuss the 'clean' correlation function K'. 
Employing the Poisson sum rule to account for the dis- 
creteness of the energy levels, one may write the disorder 
averaged density of states as 



2 



(v(E,<j>)) d = E d 2 pe^ L nmG-(p,E, ( p) = 

= T2 E f d 2 pe^e^lmG-(p,E,0), (5) 

where S&g.zxjv' ^ s an abbreviation for the sum- 
mation over all pairs k = (k x ,k y ), k x G 

2,k y £ N, L k = L(k x ,2k y ) and G~ 1 (p,E,4>) = 

(E — Tj 2 ^ — ^ (p — j^<f>e" x ) 2 ) . Substituting this expres- 
sion in the definition, Eq. (3), we obtain 



E 



■ L *l + - L *2 



2^ ^ (L k ,L k ,y/ 2 



x / e iL kl p{E F +e)-iL k2 p{E F ) _|_ c c g 



i(ki lX <j>i + k 2lX 4>2) 



(6) 



where = (L^J, p(E) = %/2mE and M = Z x M - 
(0,0). From the structure of the expression appearing 
under the angular brackets, it is clear that size and en- 
ergy averaging are equivalent: Both procedures eliminate 
contributions oscillating like exp(i ■ const. -ppL). Hence 
we obtain 



xe 



~ Lk/l [cos (k x <f> + ) + cos (k 



e— ] X 

v F , 



This expression immediately leads to 

-L k /l 



PF keM k 



Approximating the summation over the component k y by 
an integration, we finally find 



2 g-2ml/l 



(6N*(<j>))' = - PF LY, 



[cos (2m<f>) - 1] . (7) 



Technically, Eqs.(4) and (7) represent the central re- 
sult of this Letter. We next apply it to the discussion of 
the energy scales E\ , . . . , E$ . 

The first two scales E\ and E 2 are directly related 
to the ensemble averaged persistent current (of non- 
interacting electrons!) [2] 



d,L 



2 4>o 



(8) 



E\\ The curvature of the total energy at zero flux, 
i.e. the slope of the total persistent current is ob- 
tained by differentiating (8 N 2 ((f>)) ens twice with respect 
to flux (cf. Eq. (8)) [10]. As follows from Eq. (4), 

-f (i) 2 d lU=o(SN 2 (4>)) ~ ±E C , whereas the analo- 
gous contribution of (8N 2 (<j>))' is given by (cf. Eq. (7)) 



E c . We thus conclude that E\ 



is equal to E c . 

E 2 : The harmonic content of the current, is ob- 
tained by inserting the decomposition (8 N 2 (<fi)) ens = 
{8N 2 (<f>)) + {8N 2 (4>)Y into Eq. (8). Due to the large fac- 
tor pfL 3> 1 appearing in Eq. (7) the behaviour of the 
low harmonics m < L/l is dominated by the (8N 2 ((f>))'- 
contribution. For m > L/l, however, the latter is ex- 
ponentially suppressed and the harmonics take the same 
form as in the diffusive regime (cf. Eq. (1)). In particular, 
the scale governing the overall exponential suppression of 
the harmonic content is given by E c 

E3: 

dN g^,'^ one may evaluate the typical single level cur- 
rent (in the vicinity of Ep) by computing (i 2 (<f))} dL = 

{f y d i,*j Ep d ^ Kens (^ -e*,<t>u<h) ■ tws 

quantity can straightforwardly be obtained by generaliz- 
ing above derivation to the case of non-coinciding fluxes 
<f>i and fa. As a result one obtains generalizations of 
(&N 2 (4>)) respectively (8N 2 (4>))' in which the characteris- 
tic factor cos (2(f) — 1 is replaced by cos (</> + )+cos (</>_) — 2. 
Consequently we find 



A27Td</(^)) £ , d 



Employing the thermodynamic relation dI ^^ 



d<t> 



leading to E 3 = E c . 

The relevance of the Thouless energy for E\ , E 2 and 
E 3 is most easily understood from a semiclassical point 
of view: All these scales are governed by the contribution 
of high harmonics m> lto the sum Eq. (4). The behav- 
iour of any harmonic m is determined by Feynman-paths 
winding typically m-times around the ring. As soon as 
m > l/L, the motion along these paths is effectively 'dif- 
fusive' due to disorder scattering [11] and we obtain a 
corresponding result for the associated energy scales. 

Compared to these flux-dependent thermodynamic ob- 
servables, flux insensitive or transport observables change 
much more drastically in the course of a transition to the 
ballistic regime. In all cases known to us, E c loses its 
meaning and other energy scales begin to dominate. As 
an example we consider the two quantities E4 and E$. 

E4: The average two-lead conductance of a ballistic 
wire is given by the number of conducting channels lead- 



/ 



ing to E4 ~ t 

E$: It has been found in Ref. [4] that the disorder 
induced spectral correlation functions obey (orthogonal) 
Wigner dyson statistics as long as e < t _1 . We thus 
conclude E$ = t~ x (ergodicity on time scales t > r). 

In the above examples of flux dependent observables 
we have stressed analogies between diffusive and ballis- 
tic systems and emphasized the importance of disorder 
scattering. On the other hand, it is clear that all kinds 
of observables must in some or the other way approach 
their 'clean' value the more the disorder concentration 
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is reduced. In the present context it is the newly intro- 
duced function K' which guarantees the existence of a 
proper clean limit. To illustrate this, we consider a final 
example, namely the amplitude of the average persistent 
current. 
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zero flux. The linearly increasing behaviour, however, 
persists up to larger flux values <f> ~ j ^> (A/ 'E c ) 1 / 2 
leading to a maximum value J^' ~ vp/L which exceeds 
7^f x by the factor (ppL) 1 / 2 (L/l) 1 ^ » 1. For flux values 
4> > f the {I(4>)) i 

d l " curve bears resemblance with the 
sawtooth like shape obtained for clean rings [9]. 

Summarizing, we find that ballistic disordered systems 
behave to some extent 'diffusive' in their static response 
to external magnetic fluxes. In particular the scale E c , 
defined formally like in the diffusive regime, remains 
the cahracteristic energy for various observables. This 
is notwithstanding the fact that the E c in the ballistic 
regime does no longer carry the appealing significance of 
an inverse transport time through the system. When con- 
sidering physical quantities of different nature, the over- 
all picture in the ballistic regime turns out to be more 
complex than in the diffusive case. 
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FIG. 1. Average persistent current {I(<f>)) d L vs. flux 
(p F l = 1000). 

The behaviour of the persistent current in the transi- 
tion region between the diffusive and the ballistic regime 
has previously been analyzed numerically [7] but de- 
fied analytical studies so far. Earlier analytical results 
were obtained either in the strongly disordered or in the 
nearly clean limit but disagreed by a large parameter 
O(ppl), when evaluated at the boarderline L = I. Upon 
characterizing the spectrum in terms of the two func- 
tions K and K' this discrepancy disappears and one 
may smoothly interpolate between both limiting cases: 
As long as L > I, the {<57V 2 (</>))'-contribution to the 
current is exponentially suppressed (cf.fig) and we ob- 
tain the known [2,6] result for the persistent current, 

{I(<f>))^l^' in diffusive rings. Note that {li^))^^ ~ E c <f> 
increases linearly until it reaches its maximum value 
~ (AE c y/ 2 at <f> ~ (A/E c y/ 2 . Upon decreasing 
the ratio L/l, the (<5Af 2 (</>))' term gains importance: In 
the transition region L ~ I, its contribution is comparable 
to the {<5Af 2 (</>))-term and for L < I it dominates (cf.fig.). 
Like in the diffusive case, the persistent current carried 
by ballistic rings, {Ii'p))^^ se ^ s ou ^ linearly ~ E c <f> at 
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